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Abstract: We consider an extension of the i/MSM in which sterile neutrino masses orie- 
, inate from the VEV of a Higgs singlet (j) and dark matter is produced through the decays 

of (j) rather than through active-sterile neutrino oscillations. This model, which we refer 
to as the z/NMSM, can readily satisfy constraints on warm dark matter from the Lyman-a 
forest. We show that the hierarchical parameters of the z/NMSM can arise from symme- 
tries broken at or near the Planck scale for two specific models of the scalar sector: one 
in which </> stabilizes the electroweak vacuum and one in which <j) is the inflaton. These 
models have several experimental signatures that are distinct from the i/MSM, including 
additional dark radiation that is relativistic at both primordial nucleosynthesis and CMB 
decoupling and, for the former, a large invisible branching ratio of the Higgs. 
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1 Introduction 

The z^MSM [1] is an extension of the Standard Model (SM) that attempts to explain all 
observed phenomena beyond the SM using only three sterile neutrinos with Major ana 
masses below the electroweak scale. In the i/MSM, one sterile neutrino, A^i, is responsible 
for dark matter [2] while two additional sterile neutrinos, N2 and Nj,-, are responsible for 
baryon asymmetry production [3] . Moreover, the Higgs boson with a non-minimal coupling 
to gravity is responsible for inflation [4]. 

Although a detailed study of the i/MSM shows that this minimal model can explain 
most of the observed phenomena beyond the SM, there are several indications that an 
extension of the i^MSM, such as by a Higgs singlet, may be necessary: 

• Lyman-a forest bound: The Lyman-a forest [5] and other small scale structure [6, 7] 
impose severe constraints on the non-resonant production of warm dark matter in 
the z^MSM [8]. Several solutions to this problem have been proposed, including a 
resonant production of dark matter from a large lepton asymmetry [9] and a dilution 
of dark matter from a late entropy release [10-12]. Generating a sufficiently large 
lepton asymmetry requires a high level of fine-tuning or the use of a Planck-scale 
symmetry and non-renormalizable operators [13], while generating a sufficiently large 
entropy dilution requires some new physics beyond the z^MSM [10-12]. An attractive 
alternative to these scenarios uses the decays of a Higgs singlet to give a primordial 
production of dark matter [14-16] while the singlet's vacuum expectation value (VEV) 
provides an origin for the Majorana masses in the i^MSM. 

• Electroweak vacuum stability: For a Higgs mass in the range ruh ~ 125-126 GeV [17, 
18], the Higgs potential develops an instability below the Planck scale unless the 
top mass is about 2a below its central value; for its central value, an instabililty 
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develops at 10^-10^'^ GeV [19]. It has been shown that the addition of a Higgs singlet 
below the instability scale can stabilize the potential through its contribution to the 
renormalization group evolution of the Higgs quartic coupling [20, 21] or through a 
tree-level scalar threshold effect [22]. 

• Higgs inflation: There has been some discussion about unitarity violation and the 
self-consistency of Higgs inflation with a non-minimal coupling S^H^HR, where R is 
the scalar curvature and ^ ~ 10^ (see [23] and references therein). In brief, this 
model of Higgs inflation violates unitarity at the scale Aq ~ Mpi/^ when expanding 
about a small background Higgs field. Although the scale of unitarity violation is 
raised to Mpi/y/^ when expanding about the large background Higgs field during 
inflation [23], if the theory is eventually embedded into a more complete one that is 
valid up to the Planck scale then new physics is expected to appear at Mpi/S, [24]. 
This scale of new physics is below the energy scales of Higgs inflation, and therefore 
the minimality of the model is lost. Several solutions have been proposed, including 
non-renormalizable Higgs interactions that accompany the non-minimal coupling and 
restore unitarity [25] as well as strong coupling in graviton exchange processes that 
only break unitarity perturbatively [26]. However, it has not yet been shown that 
these scenarios can actually be realized [26]. Alternatively, an extension of the i/MSM 
by a Higgs singlet can "unitarize" Higgs inflation [24] ^ or provide a workable scenario 
with the singlet as the inflaton [14-16]. 

• Hierarchical parameters: The z^MSM requires a particular hierarchy in its Majorana 
masses and Yukawa couplings without an obvious origin. Motivation for such struc- 
ture can come from an extension of the underlying symmetry of the model. In the 
context of the i^MSM, models employing a U(l) flavour symmetry [27, 28], discrete 
flavour symmetries [29], and the Froggatt-Nielsen mechanism [30] have been sug- 
gested for producing the pattern of Majorana masses in the z^MSM. Such symmetries 
can provide motivation for extensions of the z/MSM involving extra singlets. 

In this paper, we consider an extension of the i/MSM by a Higgs singlet (p to address 
some of these shortcomings and argue that underlying symmetries give the required pattern 
of Majorana masses and Yukawa couplings. In particular, our starting point is a model 
in which the decays of (p allow for primordial dark matter production consistent with the 
Lyman-a forest bound. This is arguably the simplest extension of the z^MSM that provides 
an origin for the Majorana masses. We then demonstrate how symmetries broken at or near 
the Planck scale can produce the hierarchy of couplings for two versions of this model: one 
in which (f> stabilizes the electroweak vacuum [22] and one in which is the inflaton [14-16]. 

The paper is organized as follows. In Section 2, we review the constraints on the i/MSM 
and primordial dark matter production from a Higgs singlet. In Section 3, we demonstrate 
how symmetries broken at or near the Planck scale can produce the required pattern of 

^In [26], it is argued that this is not a true completion of Higgs inflation but rather an induced gravity 
inflation model added onto the SM. 
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Majorana masses and Yukawa couplings for two versions of this model. Section 4 gives the 
conclusions. 

2 The i/MSM and dark matter production from a Higgs singlet 

In this section, we first review the constraints on the i^MSM and motivate the scenario 
of dark matter production from a Higgs singlet. We then discuss the constraints on dark 
matter production from a Higgs singlet. 

2.1 The uMSM 

The Lagrangian of the z/MSM is given by 

£ = £sM + Njid^^Ni - F^iKNiH - ^NfNj + h.c, (2.1) 

where £sM is the SM Lagrangian, (/ = 1, 2, 3) are the sterile neutrinos. La {a = e, /i, r) 
are the lepton doublets, H and (p are the Higgs doublet and singlet, respectively, Faj are the 
Yukawa couplings for neutrinos, and Mjj are the Majorana masses for the sterile neutrinos. 
One of the striking features of the i^MSM is the highly constrained and hierarchical pattern 
of parameters required for successful baryogenesis and dark matter production. These 
constraints are often best expressed not in the basis Nj of (2.1) but in the basis of the 
physical mass eigenstates Nf^ with masses Mj and Yukawa couplings Fai- The two bases 
are related by the unitary transformation given in [27]. 

First, consider the constraints on and N^. The oscillations between and A'^^ 
above Tew produce a lepton asymmetry in the active neutrinos that is converted into a 
baryon asymmetry by sphalerons [3].^ and A'^^ cannot enter thermal equilibrium at 
temperatures much above Tew or else the lepton asymmetry produced in their oscillations 
is wiped out, giving the constraint [27] 

F2 < 1.2 X 10-^ (2.2) 

where Ff = (T^T)^.^. and, by convention, F2 is taken to be larger than T3 with e = T3/T2 < 
1. Similarly, masses M2, M3 <C Tew are required so that lepton number violating processes 
are negligible for T > Tew; rnasses satisfying 

M2,M3<20GeV (2.3) 

are generally considered acceptable [3, 31]. Meanwhile, effective baryon asymmetry pro- 
duction requires M2,M3 ~ M to be highly degenerate with a mass difference AM = 
M3 — M2 <C M [3]. The baryon asymmetry produced can be expressed as a function of 
F2,€,M,AM, and the neutrino hierarchy. Since active neutrino masses are generated via 
the seesaw mechanism, one of F2, e, and M (typically F2) can be expressed in terms of the 
others with the relation [27] 

Arriatm ^ -, (2.4) 

2M 

^Tew ==: 140 GeV for a Higgs mass mn ~ 125 GeV [31]. 
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where Amatm — 0.05 eV, v = 246 GeV, and k = 1(2) for the inverted (normal) hierarchy. 
Analytic expressions for the baryon asymmetry are given in [3, 32] while a numerical study 
has been carried out in [31]. The allowed range of each parameter individually is [31] 

M > 140 MeV, (2.5) 
10"^ eV < AM < MeV, (2.6) 
10-4 < e < 1, (2.7) 



though the combination must produce the observed baryon asymmetry hb/s ~ (8.4-8.9) x 
10-11 [33]. Note that the lower bound (2.5) comes from demanding that and A'^^ 
decay before Big Bang nucleosynthesis (BBN) [34, 35]'^ and that a significant amount of 
parameter space for M < 500 MeV is ruled out by the CERN PS191 experiment and other 
accelerator bounds [31, 35]. For the parameter space allowed by (2.5)-(2.7), a lower bound 
on F2 is approximately 

F2 > 3 X 10-^ (2.8) 

Now consider the constraints on the dark matter candidate N^. Mixing with active 
neutrinos below Tew allows the 1-loop decay A™ — )• u'^j with width [1, 37] 



9aG? 



5 



(2.9) 



where Of = v^F^/ (2Mf) and is evaluated with [27] 

Fal ~ Tal + ^F^2 + ^Fa3. (2.10) 

The second and third terms on the right hand side of (2.10) are contributions to F^i 
induced by the mixing of Ai with A2 and N3 to form the mass eigenstate A™. Direct 
searches for the X/7-ray line corresponding to this decay provide the strongest limits on 
9i (as a function of Mi) for the mass range relevant to the z/MSM. A summary of these 
limits is given in [8]. In general, 

^?<3xl0-5(^y (2.11) 

must be satisfied for 0.5 keV < Mi < 14 MeV, though the constraint is typically 100 
times stronger than (2.11) for masses outside the 12-40 keV range [8]. For A{" produced 
entirely from active-sterile neutrino oscillations, Mi can be bounded above by combining 
the X-ray constraints with the requirement of sufficient dark matter production (oc 9l). 
The bound obtained depends on the lepton asymmetry at the time of A^J" production: a 
negligible lepton asymmetry is called the non-resonant production (NRP) scenario while a 



^Recent work [36] suggests this bound can been strengthened to M > 1.4 GeV in tlie RP scenario, 
discussed later. 
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large lepton asymmetry is called the resonant production (RP) scenario. The bounds for 
these two scenarios are [8, 9, 11] 



Mf^P < 2.2 keV, < 40 keV. (2.12) 

Meanwhile, Mi can be bounded below by phase-space density arguments for dwarf spheroidal 
galaxies [38], the Lyman-a forest data [5, 39], studies of gravitationally lensed QSOs [6], 
and N-body simulations of the Milky Way [7]. The bounds from the Lyman-a forest data 
and N-body simulations of the Milky Way are the strongest and give^ 

^NRP > ^3 j^gy^ ^RP > 2 ]^gv. (2.13) 

Combining (2.12) and (2.13) rules out the simpler NRP scenario, even with a possibly large 
entropy dilution from the decays of A^2" ™d [11]. The RP scenario is still allowed for a 
range of Mi; it requires an even larger degeneracy than (2.6), on the order AM < 10~^ eV, 
to produce the required lepton asymmetry for enhanced dark matter production [13]. This 
level of degeneracy is unstable in the presence of radiative corrections and must be achieved 
with either fine-tuning or an extension of the model by a Planck-scale symmetry and non- 
renormalizable operators [13]. 

2.2 Dark matter production from a Higgs singlet 

An alternative dark matter production scenario that is capable of satisfying the Lyman-a 
forest bound uses a real Higgs singlet (j) and its decays to N^^ [14]. This scenario is arguably 
simpler than the RP scenario and has the advantage that Majorana masses originate from 
the VEV of (p rather than as bare mass terms. This extension of the i^MSM, which we will 
call the neutrino Next-to- Minimal Standard Model (z^NMSM), is the basis of this paper. 
In the z^NMSM, the Majorana mass term in the Lagrangian (2.1) is modified to 

A£ = -^,/.iVfiVj, (2.14) 

where M/j = Xjj (cp) once (p acquires a VEV. In the mass basis Nf^, A/j — )• A/ where 
M/ = Xj {(p). The mixing angle 6f is assumed small enough that dark matter production 
from active-sterile neutrino oscillations is negligible and the X/7-ray constraint (2.11) is 
satisfied.^ Assuming no miraculous cancellations of terms in (2.10), this requires 

F.i,^F.„^F„3<10-^ (2.15) 

With (2.15), one can show that the induced contributions to Mi from M12 and M13 are 
small [27] and hence Ai ~ An. Dark matter production then proceeds via the decays 



■'These are the (Bayesian) 2a bounds. Although [7] quotes a stronger bound for Aff , it is based on a 
simple mass rescaling argument that is shown to be insufficient for a more rigorous analysis of the Lyman-a 
forest bound in the RP scenario [39]. 

^Since the NRP and RP bounds (2.12) no longer apply, Mi may exceed the range in which (2.11) is 
valid. In this case, 7-ray constraints from EGRET [40] and FERMI [41] give r > 10^^ s, or equivalently 
0l<2x 10-2° (MeV/Ml)^ for masses up to 30 TeV. 
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_^ jvf 7V[" with the partial width [14] 

\2 ^2 

r = ~ ^ms, (2.16) 

levr ^ 167r ^ ' 

where > 2Mi is the mass of the physical mass eigenstate c^™ .^ This production depends 
on the thermal history of cj)^, specifically the ratio of its mass to its freeze-out temperature, 
Tf = m^/Tf [42]. For the case that is in thermal equilibrium down to temperatures 
T <^ rriffi (i.e. rj ^ 1), the dark matter abundance is given by [14] 

'J^IS-^ " , (2.17) 

where / (m^) ~ (10.75/c/* (m<^/3))^/^ and 1 < S* < 2 is a factor that accounts for entropy 
production from the decays of and iV™ after A^"^ is produced.'^ Using Mi ~ An {(j)) in 
(2.17), the appropriate dark matter abundance $7^?™ ~ 0.23 is generated when 



For the case that (f/^ is a thermal relic decaying out of equilibrium (i.e. rj <^ 1), the dark 
matter abundance is given by [42] 

0.3 / Ml \ / 10.75 \ f B \ , , 

where B = F/F^"* is the branching ratio of (j)"" N^N^.^ Analytic expressions relevant 
to the intermediate case rj ~ 1 can be found in [42], and the result is a combination of 
(2.17) and (2.19). 

The Lyman-a forest bound for this dark matter production mechanism can be esti- 
mated by rescaling the NRP bound, giving [11, 42] 

1 /3 

M^'"' > 10 f keV, (2.20) 

where Tprod is the temperature at which N^^ is produced. Further constraints come from the 
requirement that the interactions (/)™ o N^N^ and •(->■ N^N^ (and any interactions 
SM ^ N^N^ and SM o A^f A^g™ mediated by (/>™) do not bring and into thermal 
equilibrium at the characteristic temperature of leptogenesis [43] 

Tl ~ ^ j , (2.21) 



®We have assumed a small mixing angle Oh<i> between the Higgs boson h and 4> so that 4>"^ ~ (f>, h"^ ~ h, 
and the decays h"^ — >• N^^N^ are negUgible compared to (ji™ — ^ N^N^ [42]. This is a good approximation 
for both models considered in Section 3. 

''Since A'^^ production peaks at Tprod = 71^^/2.3, (2.17) is a good approximation for r/ > 3 [42]. 

*As in [42], we neglect any (jxj) annihilations that could reduce (2.19) by up to a factor of 2. 
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where Mq ~ 7 x 10^^ GeV, and spoil baryogenesis.^ Moreover, the addition of (j) must 
not open an invisible branching ratio of the Higgs greater than 65% at 2a [44]. These 
constraints are discussed further in Section 3 for specific models of the scalar sector. 

Although we have assumed that (p is real in the discussion above, the possibility of a 
complex (p is worth considering since it allows for a richer symmetry structure to give the 
masses and couplings of the i/NMSM. (We parametrize (f> = {p + ix) / for a complex (j) 
but use m<^ and 0™ instead of nip and /O™ to maintain consistency with the notation for a 
real (j).) Typically, cj) is assumed real to avoid a massless Goldstone boson x and hence the 
unsuitably fast decay channel N^^' — t- z^™x for dark matter [14, 42, 45]. However, there are 
several ways to avoid this decay for a complex (p. First, if (j) is charged under a discrete 
symmetry and terms of the form (p"^ + (<^^)" are allowed, these terms give x a mass and 
can kinematically forbid the decay N^^ — ?• i^"^X- If analogous decays N^^^N'!^ — t- v'^x 
are still allowed, they can relax the constraint (2.5) to M > few MeV [27]. Alternatively, if 
(p is charged under a global continuous symmetry, ^'^ the mixing of A'^i with other neutrino 
species can be sufficiently reduced (or even forbidden) by a symmetry to suppress the decays 
A^"^ — )• to a massless x- This case is particularly interesting since x can contribute 
to the effective number of neutrino species and give a value of N^s above the SM (or 
z/MSM) prediction, as recent measurements prefer (see [46] and references therein). The 
contribution of x to Aeg depends on the freeze-out temperature Tf. it can be as large 
as AAcff ~ 1 for a thermal distribution of x or much smaller if x decouples early. The 
Planck experiment and other future cosmic microwave background (CMB) experiments 
will therefore be able to constrain these models with a complex (p [48]. 

3 Symmetries and the j/NMSM 

The i^NMSM, like the z/MSM, requires parameters that are constrained to be hierarchi- 
cally small. In the context of the i^MSM, flavour symmetries and the Froggatt-Nielsen 
mechanism have been explored for producing the required pattern of Majorana masses and 
Yukawa couplings [27-30]. Following this approach, we would like to demonstrate how the 
parameters of the z^NMSM can arise from symmetries broken at or near the Planck scale. 
Since the values of some parameters (e.g. An in (2.18)) depend on an unspecified scalar 
sector, we first keep the discussion general and then consider two specific models of the 
scalar sector: one in which (p stabilizes the electroweak vacuum [22] and one in which (p 
is the inflaton [14-16]. These models of the scalar sector, though motivated as minimal 
solutions to shortcomings of the i/MSM, are only meant to provide definite examples for 
the symmetries used in the flavour sector; other models may certainly be considered. We 
do not provide an explanation for the values of parameters in the scalar sector or the 
associated hierarchy problem since little is known about their origin. 

®If these interactions bring and into thermal equilibrium below Tz,, the asymmetry in the sterile 
neutrinos will be destroyed but the asymmetry in the active neutrinos will remain. 

^"We assume a global symmetry to avoid introducing a new low-energy gauge sector. 

^^The real component of (p can also contribute to A^ofi during BBN if < 10 MeV [47]. For the models 
of the scalar sector considered in Section 3, however, 3> 10 MeV. 
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Ni N2 N3 H <p 



U(l) 
Z3 



-1 -1 -1-10 2 



1 



1 



1 



1 







1 



Table 1. Examples of an anomaly-free global U(l) and Z3 symmetry that can give the Lagrangian 
structure (3.1). The SM fermions not shown have B — L charges under the U(l) while only Ea (the 
right-handed charged Icptons) have charges, 1, under the Z3. 

First consider how the structure of the i/NMSM Lagrangian, 



can arise from an underlying symmetry without regard to the size of the couplings Faj and 
A/ J. There are several ways this structure can arise: 

• Conformal symmetry /scale invariance: The structure (3.1), which has only terms 
with dimensionless couplings, can arise from models with a classical conformal sym- 
metry [49] or hidden scale invariance [50]. These models have been motivated as a 
solution to the hierarchy problem: the conformal symmetry forbids tree-level scalar 
mass terms while radiative breaking of this symmetry by the conformal anomaly is 
responsible for electroweak symmetry breaking and, in [50], a hierarchy between the 
electroweak and Planck scales from a choice of large scale /. Unfortunately, existing 
models of this type are not fully realistic. 

• Global U(l) symmetry: For a complex 0, the structure (3.1) can arise from a global 
U(l) symmetry under which (j) is charged (we use a global symmetry to avoid in- 
troducing a new low-energy gauge sector). Such symmetries can arise from string 
theory as the remnant of a non-linear ly realized U(l) gauge symmetry in which the 
gauge boson acquires a large (string scale) mass through its coupling to a Stueckel- 
berg field [51]. For a consistent model, the underlying U(l) gauge symmetry must 
be anomaly-free or Green-Schwarz anomalous [52, 53]. An anomaly-free example in 
which matter fields have U(1)b-l charges is given in Table 1. 

• Discrete Zjsi symmetry: A discrete Z]\j symmetry can also give the structure (3.1). 
Such symmetries can arise from the spontaneous breaking of a continuous symmetry 
at a high scale [54] or from coupling selection rules on heterotic orbifolds (see [55] and 
references therein). Note that it is often easier to satisfy the anomaly cancellation 
conditions for Z^r symmetries [55, 56] than those for U(l) symmetries: an anomaly- 
free Z3 example is given in Table However, the spontaneous breaking of discrete 
symmetries can produce domain walls [59] and care must be taken to avoid these, 
such as by breaking the discrete symmetry before inflation [60]. 



^^The mixed .Zjv-U(l)i'-U(l)y anomaly does not pose a meaningful constraint since the hypercharge 
normalization is not fixed [57, 58]. 
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^ 




(0 F,2 Fe3\ 


A/j = 





A22 A23 


, Fal = 


F,,2 ^^3 






A23 ^33^ 




[0 Fr2 Fr3 J 



Although either a global U(l) or discrete Z^r symmetry can give the desired Lagrangian 
structure (3.1), we use a global U(l) to avoid introducing domain walls. 

Now consider the hierarchy of Majorana masses and Yukawa couplings in the i^NMSM. 
To explain the small Yukawa couplings Fai ^ 10^^^ and, for a complex <p, to prevent the 
fast dark matter decay channel — )• I'^Xy we introduce a Z2 symmetry under which 
only A^i is charged (see Table 2).""^^ This symmetry allows only the couplings 



(3.2) 



and hence forbids mixing of A^i with the other neutrinos, making N^^' completely stable 
(^1 = 0) and one active neutrino exactly massless. The required pattern of Majorana 
masses and Yukawa couplings can then be produced if there are strong hierarchies in the 
remaining A/j and Fai, specifically if 

All ~ TTT' A23 ~ 7TV' ™ax|A22, A33j ~ , 

(W (W {9) (3.3) 

Fa2 ~ F2, Fa3 ~ F3. 

We consider two possibilities for generating these hierarchies: 

• Frog gatt- Niels en mechanism: The Froggatt-Nielsen mechanism [61] is a well-known 
method of generating hierarchical parameters. In brief, a new U(1)fn gauge symme- 
try that is spontaneously broken by a flavon field "i? at a very high scale is introduced. 
Fields of the i/NMSM are charged under this U(1)fn so that (or -d^) must couple 
to the terms in (3.1) with various powers to form gauge singlets. After the U(1)fn 
is spontaneously broken, these non-renormalizable terms are suppressed by powers 
of r/ = (i?) /Mpi, where r/ is a free parameter (though typically assumed to be on 
the order of the Cabibbo angle [30, 62]). Of course, multiple flavon fields i^j with 
various rji = {-di) /Mpi may be used, as well as a discrete Zjy symmetry in place of 
the U(1)fn. 

• Non-perturbative symmetry breaking: Another possibility for generating hierarchical 
parameters comes from non-perturbative symmetry breaking in string theory. In [63] , 
for example, it is shown that heterotic string compactifications on Calabi-Yau man- 
ifolds can give models with the SM gauge group and additional U(l) symmetries. 
These additional symmetries can play a role analogous to that of the U(1)fn: if 
the z^NMSM fields are charged under these symmetries, the terms in (3.1) may re- 
quire couplings to various powers of "di = e~'^'^^^^^^ to form gauge singlets, where 

= f + 2ix^ are Kahler moduli with axionic components (not to be confused 
with the Goldstone boson x) that transform non-linearly under the U(l). After these 
symmetries are spontaneously broken by (t*) ^ Afpi [63, 64], the terms in (3.1) are 



^^The anomaly cancellation conditions for this Z2 are trivially satisfied. Therefore this symmetry is exact 
at the quantum level. 
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suppressed by powers of r/j = e ^* )/*'^pi. Again, discrete symmetries may be used 
in place of the U(l) symmetries. 

Although either mechanism may be used to generate the hierarchical parameters (3.3) for 
the same charge assignment, the non-perturbative symmetry breaking does not require 
additional symmetry breaking or scalar particles below the Planck scale and hence adheres 
closer to the "minimal" philosophy of the j^MSM. From a more practical point of view, the 
non-perturbative symmetry breaking may be necessary if a symmetry, used in place of 
U(1)fn5 would otherwise be broken after inflation and produce domain walls [60]. 

To fix the absolute scale of the couplings A/j, the values of and {(j)) must be fixed 
(see (2.18) and (3.3)) by some model of the scalar sector. We now consider two models 
of the scalar sector that are motivated as solutions to the shortcomings of the z^MSM, as 
discussed in Section 1. 



3.1 Stabilization of the electroweak vacuum 

A model of the scalar sector in which a Higgs singlet can stabilize the electroweak vacuum 
through a scalar threshold effect is given in [22]. This model considers a complex (j) and 
scalar potential of the form 



V 



V = Xh [Hm - - 
+ 2Xh^ ( WH 



w 

T 

2 

W 

T 



(3.4) 



which is the most general renormalizable potential that respects a global abelian symmetry 
under which only (j) is charged. Values of A/,,, A^ > and < Xh^cj) are assumed so that 
the minimum of this potential is given by 



w 

T 



(3.5) 



where v = 246 GeV. The mass matrix for the real components of H and (p is then 



( Xhv"^ Xh^vw\ 

yXh^vw x^w"^ J ' 



(3.6) 



while the imaginary component of (f) (i.e. x) remains massless. In contrast to other models 
that use a Higgs singlet to stabilize the electroweak vacuum (e.g. [20, 21]), this model 
assumes w ^ v. The two eigenstates of (3.6) then have masses 



rrij 



2v' 



2w^ 



X2 



+ 



X2 



X6 + ^^ 



(3.7) 
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with a mixing angle 9^^ ~ v/w. Integrating out the heavier state for scales below gives 
the effective potential 

V,s = \{h^H-'^^ , A^A^-^, (3.8) 

where the matching condition for the Higgs quartic coupling gives a tree-level shift 5\ = 
Xj^^/Xff) from A just below to \h just above m^. Provided is below the instability 
scale A ~ lO^-lO^o GeV [19], a value of SX ~ 0.01 can push the instability beyond the 
Planck scale. 

Due to the massless Goldstone boson X) the value of A/j,^ is constrained by limits on 
the invisible branching ratio of the Higgs. For nifi ~ 125 GeV, the total SM decay width 
of the Higgs is [65] 

TsM = 4.07 MeV, (3.9) 



while the invisible decay width for /i™ — ?• xx is [20] 



(3.10) 



Allowing an invisible branching ratio of up to 65% [44] gives the constraint 

Ah<^ (m^,) < 0.02. (3.11) 

A value of SX that stabilizes the electroweak vacuum and is consistent (3.11) can then be 
obtained for A^ < 0.04 (the running of Xhcf, and A,^ is small for these values). We illustrate 
this by constructing a model with A/i^,A^ ~ 0.01. For these values, there is an invisible 
branching ratio of the Higgs of about 30%. Moreover, one can show that x remains in 
thermal equilibrium down to temperatures just below m^. The model therefore has a 
AA'cff ~ 4/7 contribution to the effective number of neutrino species from x and hence a 
total value of N^s — 3.6. This value can be tested by the Planck experiment and other 
future CMB experiments [48]. 

Now consider the flavour sector of the i^NMSM for this model of the scalar sector. For 
Xhtp ~ 0.01, the interactions H^H -H- (^"^(j)^ keep (jf^ (the real component of (p) in thermal 
equilibrium down to temperatures T <^ for any mass < A. We are therefore in the 
dark matter production case rj ^ 1. For A^ ~ 0.01, the ratio m^/ {(f)) is fixed by (3.7) and 
(2.18) gives a value of 

All ~ 1 X 10"^ (3.12) 

to produce the correct dark matter abundance. The Lyman-a forest bound (2.20) is 
therefore satisfied for a choice {(f)) > 500 GeV. Taking (i;^) ~ 10^ GeV, a pattern of masses 
M/j and couplings F^i that gives the correct dark matter abundance and baryon asym- 
metry can be achieved with two fields '!?i,t?2; the values rji ~ 10~^,r72 ~ 10^^, and the 
charge assignments given in Table 2. We stress that this is the simplest example of an 



^^Here we have used 5 ~ 1 due to the fast decays 7V™3 — >■ 8.nd taken > Tew for / (m^). Also, 
All must be a factor of larger than in (2.18) for a complex since only the real component of 4> can 
decay to iVf. 



- 11 - 





Ni 


N2 


A^3 


La 


H 




^1 


^2 


U(l) 


-1 


-1 


-1 


-1 





2 








u(iy 

















-1 


1 










1 


-1 














1 


Z2 


1 
























Table 2. Charge assignments for the stabihzation of the electroweak vacuum scenario. The SM 
fermions not shown have B — L charges under the gfobal U(l) and are uncharged under the global 
U(l)', Z3, and Z2. 
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anomaly- free model we could find, though other charge assignments are possible. 

In this model, we use a global U(l) symmetry to give the structure (3.1) and the 
potential (3.4), as well as a global U(l)' and discrete Z3 and Z2 symmetries to give the 
textures 

/r?i \ (^mm\ 

(3.13) 

(Note that we implicitly assume the terms 'dl'diH^H and 'd^'dicl)'^ (p are sufficiently small to 
preserve (3.4).) The parameters of the z^NMSM are then 

Ml ~ 1 GeV, M ~ 1 GeV, AM ~ 100 eV, 
F„2 ~ 1 X 10-^ F«3 ~ 1 X 10-^ 

up to O (1) constants. Active neutrino mixing is anarchical (up to charged lepton correc- 
tions) while the charged lepton and quark Yukawa couplings remain unsuppressed. There- 
fore additional flavour symmetries with anomaly cancellation techniques along the lines 
of [30, 66] must be used to produce hierarchies in the charged lepton and quark sectors. As a 
consistency check on this model, we must verify that A*"!" and A3" are out of thermal equilib- 
rium at the leptogenesis temperature Tj, ~ 3 x 10^ GeV. Since ~ 2 x lO'^ GeV Tl, cj/^ 
has decayed away by Tl and we need only consider the scattering processes H'^ H o N^N^ 
and o N^N^ mediated by cj)"" and x-^^ For Xh^ ~ 0.01, these processes are out of 

equilibrium at for A2,A3 ~ A23 ^ 10^^, which is certainly satisfied by (3.13). 

This model demonstrates that it is possible to use symmetries broken at or near the 
Planck scale to obtain the hierarchical pattern of Majorana masses and Yukawa couplings 
required for successful baryogenesis and dark matter production. The model obeys all 
phenomenological constraints and allows for the possibility of Higgs inflation by ensuring 
that the Higgs potential does not develop a second minimum before the Planck scale. Of 



^^The Z3 and Z2 symmetries could be combined in a single Ze, if desired. 

'^^Note that a relic CP-even distribution of and N^^ from the decays of (/>™ does not affect leptogenesis. 

^^If the Froggatt-Nielsen mechanism is used, the spontaneous breaking of the Z3 symmetry below the 
Planck scale may or may not lead to domain walls (depending on the underlying potential) [60]. If necessary, 
a simple way to avoid these domain walls is to take 772 smaller (i.e. break the discrete symmetry at a higher 
scale) and require higher powers of to couple to N2 and N3. For example, one could take 772 — 5 x 10^'^ 
and replace the Z3 with a Zg under which (A^2,A^3,i?2) have charges (3,-3,1). 
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course, the symmetries used do not address the hierarchy problem associated with radiative 
corrections to the Higgs sector. To do so would involve implementing a supersymmetric 
version of the theory, which departs from the underlying philosophy of the i^MSM, or 
implementing a conformal symmetry solution, which requires an understanding of how to 
include gravity in such a theory. This is something we cannot do at present. 

3.2 (p Inflation 

A model of the scalar sector in which (j) can be a light inflaton (m^ < ruh) and thus provide 
an alternative to Higgs inflation is given in [14] and developed further in [15, 16]. The 
scalar potential of this model is 

V {H,ct)) = A (^H^H - jct^Uf + J {ckUy - \rn^(l>U, (3.15) 

where m <^ ^/]3M■p\ so that chaotic inflation proceeds via the quartic term and, in contrast 
to [14-16], we have taken cj) to be complex so that the Lagrangian structure (3.1) can arise 
from a global U(l) symmetry. The potential (3.15) is then the most general renormalizable 
potential that respects a global U(l) symmetry under which only (p is charged, assuming 
the bare mass term for the Higgs is negligible. Successful chaotic inflation^^ requires /? ~ 
1.5 X 10"^'^ to give the correct amplitude of adiabatic scalar perturbations and a < 10"'^, 
A/j < 1.5 X 10"^ to ensure that the flatness of the potential is not spoiled by radiative 
corrections from the loops of SM particles and sterile neutrinos [16]. Achieving a sufficiently 
high reheating temperature Tj. > for baryogenesis requires a >7x 10"^*^ [15]. Moreover, 
a value of A ~ 0.13 is required for rrih ~ 125 GeV. For these parameters, expanding the 
potential (3.15) about its minimum gives the relations 




(3.16) 



where v = 246 GeV. Note that we do not provide an explanation for the small values of a 
and /3 in the scalar potential; we simply take their values to be within the range allowed 
by successful inflation. Also note that, for a < 10"^, the invisible branching ratio of the 
Higgs is negligible. 

Now consider the flavour sector of the z^NMSM for this model of the scalar sector. 
As in [14], we assume an inflaton mass 300 MeV < < nifi so that the mixing angle 
9fi^ is large enough to keep (p"^ in thermal equilibrium down to temperatures T <^ 
via the interactions 0™ -H- e~e'^,(f)"^ -H- etc.. We are therefore in the dark matter 

production case rj » 1. The ratio m^/ {(p) = y/]3 is fixed by (3.16) and (2.18) gives a value 
of 

All ~ 3 X 10"" (3.17) 

^*For a real 0, (3.15) was originally presented as the most general scale-invariant potential in which the 
scale invariance is explicitly broken by a mass term for [14]. 

^^As mentioned in [14], chaotic inflation with a quartic potential is disfavoured by WMAP data [67]. 
However, a non- minimal coupling to gravity can help bring this model in line with the data [68]. 
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Table 3. Charge assignments for the (j> inflation scenario. The SM fermions not shown have B — L 
charges under the global U(l) and are uncharged under the global U(l)', Z4 and Z2. 



to produce correct dark matter abundance. The absolute scale of ((/>), however, is not 
fixed. There is a relatively narrow window 7 x 10^ GeV < ((/>) ^ 5 x 10^ GeV that is 
consistent with the constraints on a, the assumption > 300 MeV, and the Lyman-a 
forest bound. Taking {cj)) ~ 10^ GeV, a pattern of masses M/j and couplings Faj that 
gives the correct dark matter abundance and baryon asymmetry can be achieved with two 
fields ??i,'i?2) the values ryi ~ 2 x 10"'^, 172 — 5 x 10^^, and the charge assignments given in 
Table 3. Again, this is the simplest anomaly-free model we could find, though other charge 
assignments are possible. 

In this model, we use global U(l) and U(l)' symmetries to give the structure (3.1) and 
the potential (3.15) (the U(l)' forbids the terms (j)'^{'d\)'^ + {(ji^^^d^ and hence prevents a 
discrete Z3 from forming once 7?i acquires a VEV), as well as discrete Z4 and Z2 symmetries 
to give the textures 







r]fii 



rjiri2 r/ir/2 
rjiri2 r/ir/2 
\0 rjiri2 r/i?72^ 



(Again, we implicitly assume the couplings for 'd^'diH^H and 
to preserve (3.15).) The parameters of the z^NMSN are then 

Ml ~ 20 keV, M ~ 4 GeV, AM ~ 10 eV, 



F„2 ~ 1 X 10" 



'a3 



1 X 10" 



(3.18) 



are sufficiently small 



(3.19) 



up to O (1) constants. As before, active neutrino mixing is anarchical (up to charged lepton 
corrections) and additional flavour symmetries must be used to produce hierarchies in the 
charged lepton and quark sectors. We also have the parameters 



a 



4 X 10"^ ~ 300 MeV, Oh^c^lx 10" 



(3.20) 



For these values, it can be shown that x remains in thermal equilibrium roughly while 
does (to temperatures below m^) via the interactions (p"^ -H- XX- The massless x therefore 



^"We have anticipated 5 ~ 1 and ~ 300 MeV in obtaining (3.17), though these parameters only have 
an O (1) effect on An. 

^^For ruh — 125 GeV, the constraints on a can be strengtliened by fimits on axion searches in the CHARM 
experiment [16]. The mass range allowed by this experiment, 270 MeV < < 1.8 GeV, corresponds to 
2 X 10"^" < a < 8 X 10"^ 
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contributes AA'gff ~ 4/7 to the effective number of neutrino species. As a consistency 
check on this model, we must verify that NI^ and are out of thermal equilibrium at 
the leptogenesis temperature T/^ ~ 2 x 10'^ GeV. Since < 2M, the processes (p"^ — >• 
N^N^ and cj)^ — )• N^N^ are kinematically forbidden and the dominant processes are 
H^H o N^N^"- and H^H o N^Nr/'. These are out of equilibrium at Tl for A23 < 0.01, 
which is satisfied by (3.18). One can also verify that the reheating temperature can be 
as large as T^. ~ 5 x 10'^ GeV for a ~ 4 x 10~^ [15], which is above the leptogenesis 
temperature. 

This model demonstrates that, for a scenario in which <^ is a light inflaton, it is again 
possible to use symmetries broken at or near the Planck scale to obtain the pattern of Ma- 
jorana masses and Yukawa couplings required for successful baryogenesis and dark matter 
production. This model obeys all phenomenological constraints and provides an alternative 
to the Higgs inflation of the z^MSM, but it does not improve the stability of the electroweak 
vacuum. 

4 Conclusion 

The i/MSM is an extension of the SM that attempts to explain neutrino oscillations, dark 
matter, the baryon asymmetry of the universe, and inflation using only three sterile neu- 
trinos with masses below the electroweak scale. Despite the phenomenological successes 
of the z^MSM, an extension of the model may be necessary to accommodate the Lyman-a 
forest bound, stabilize the electroweak vacuum, allow for Higgs inflation, and provide an 
explanation for its hierarchical pattern of masses and couplings. In this paper, we have 
studied how an extension of the j^MSM by a Higgs singlet 4> can address these issues. 

Our starting point has been a model in which the decays of (p give primordial dark 
matter production consistent with the Lyman-a forest bound and in which the VEV of (p 
produces the Major ana masses of the sterile neutrinos. From this next-to-minimal model, 
or i/NMSM, we have considered two versions of the scalar sector that can address other 
shortcomings of the z^MSM: one in which a complex (p stabilizes the electroweak vacuum 
via a scalar threshold effect, thereby allowing for the possibility of Higgs inflation; and 
one in which a complex (p provides an alternative to Higgs inflation but does not improve 
the stability of the electroweak vacuum. For both of these models, we have demonstrated 
how symmetries broken at or near the Planck scale can produce the required hierarchical 
pattern of Major ana masses and Yukawa couplings. The former model uses a U(l)^ x Z3 x Z2 
symmetry while the latter uses a U(l)^ x Z4 x Z2 symmetry. We have not, however, provided 
an explanation for the parameters of the scalar sector or addressed the hierarchy problem 
associated with radiative corrections to the Higgs sector. 

The models presented in this paper satisfy all phenomenological constraints and make 
several experimental predictions that are distinct from the i/MSM. These predictions in- 
clude completely stable N^^ dark matter (hence no visible X/7-ray signals from its decays) 
and anarchical active neutrino mixing angles (up to charged lepton corrections) with one 
active neutrino exactly massless. Moreover, both models have N^fi — 3.6 for the effective 
number of neutrino species while the former model has an invisible branching ratio of the 
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Higgs of about 30%. It will therefore be possible to test these models with the Planck 
experiment and the LHC in the near future. 
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